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Abstract 

We have included the effective description of squark interactions with chargi- 
nos/neutralinos in the MadGraph MSSM model. This effective description includes 
the effective Yukawa couplings, and another logarithmic term which encodes the 
supersymmetry-breaking. We provide an extensive test of our implementation ana- 
lyzing the results of the partial decay widths of squarks into charginos and neutralinos 
obtained by using FeynArts/FormCalc programs and the new model file in Mad- 
Graph. We present results for the cross-section of top-squark production decaying 
into charginos and neutralinos. 



1 Introduction 



The Standard Model (SM) of the strong and electroweak interactions is the present pa- 
radigm of particle physics. Its validity has been tested to a level better than one per 
mile at particle accelerators pQ. Nevertheless, there are arguments against the SM be- 
ing the fundamental model of particle interactions [2], giving rise to the investigation of 
competing alternative or extended models, which can be tested at high-energy colliders, 
such as the Large Hadron Collider (LHC) [2111], or a 500 — 1000 GeV, e + e~ International 
Linear Collider (ILC) [SJE]- One of the most promising possibilities for physics beyond the 
Standard Model (SM) is Supersymmetry (SUSY) [7HTU]. which leads to a renormalizable 
field theory with precisely calculable predictions to be tested in present and future exper- 
iments. The simplest model of this kind is the Minimal Supersymmetric Standard Model 
(MSSM). Among the most important phenomenological consequences of SUSY models, 
is the prediction of new particles, the SUSY partners of SM particles, sometimes called 
sparticles. There is much excitement for the possibility of discovering these new particles 
at LHC lUH]. The LHC collaborations are already performing searches on these particles, 
and excluding portions of the SUSY parameter space, see e.g. [T2H27] . Recently, the CMS 
and ATLAS collaborations have reported a 5 standard deviations signal on a new boson 
particle, at a mass m ~ 125 GeV, which is compatible with the interpretation of a SM 
Higgs boson [2"SH3"T] . The CDF and DO collaborations at the Tevatron also found (less sig- 
nificant) signals compatible with this interpretation [32, 33] . This signal is also compatible 
with the lightest neutral Higgs boson of the MSSM (see e.g. [31]). Precision measurements 
and precision computations are both mandatory nowadays. Under the huge avalanche of 
new data at present, an accurate prediction of sparticles couplings to other particles and 
their production cross-section is needed. In this work we focus on the properties of the 
squarks - the SUSY partners of SM quarks. In particular, we concentrate on the squark 
decay channels involving charginos and neutralinos - the fermionic SUSY partners of the 
electroweak gauge and Higgs bosons. 

Once produced, squarks will decay in a way dependent on the model parameters (see 
e.g. [35] )• If gluinos (the fermionic SUSY partners of the SM gluons) are light enough, 
squarks will mainly decay into gluinos and quarks (q — > qg) [36|I3"T]. which proceeds trough 
a coupling constant of strong strength. If the mass difference among different squarks is 
large enough, some squarks can decay via a bosonic channel into an electroweak gauge 
boson and another squark (q a — > q' b (Z, H /± )), and if Higgs bosons are light enough, also 
the scalar decay channels are available (q a — > q' b (h°, H°, A , H^)) [3HrlH], which can be 
dominant for third generation squarks due to the large Yukawa couplings. Otherwise, the 
main decay channels of squarks are the fermionic ones: chargino/neutralino and a quark 
(q —> q'x) [12H11]. Some of those channels are expected to be always open, given the 
large mass difference between quarks and squarks, and that the charginos/neutralinos are 
expected to be lighter than most of squarks in the majority of SUSY-breaking models. In 
the few cases in which these channels are closed, the squarks will decay through flavour 
changing neutral channels [13H1E], or through three- or four-body decay channels involving 
a non- resonant SUSY particle [IHHSl]- 
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Here we will concentrate on the squark decay channels involving charginos and neutra- 
linos. Their partial decay widths were computed some time ago, including the radiative 
corrections due to the strong (QCD) [55H57]. and the electroweak (EW) [321 1331 1551 155] 
sectors of the theory. These radiative corrections are large in certain regions of the param- 
eter space [12], and their complicated expressions are not suitable for their introduction 
in the Monte-Carlo programs used in experimental analysis. Recently Ref. [HJ presented 
an improved description of squark/chargino/neutralino couplings, more simple to write 
and to introduce in computer codes. This computation combines the effective descrip- 
tion (which includes higher order terms) with the complete one-loop description (which 
includes all kinetic and mass-effects factors) and defines a new effective coupling. It in- 
cludes a non-decoupling logarithmic gluino mass term, which implies a deviation of the 
higgsino/gaugino and Higgs/gauge couplings equality predicted by exact SUSY. This de- 
viation is important and has to be taken into account in the experimental measurement 
of SUSY relations. Ref. [H] showed that the effective description approximates the im- 
proved description within a 10% precision, except in special uninteresting corners of the 
parameter space, where the corresponding branching ratios are practically zero. Ref. [HJ 
applies the description only to squark decays. The present work expands the results of 
Ref. [HJ by applying those results to the production cross-section of squarks at the LHC. 
To that end we have implemented this effective description in MadGraph's [601462] MSSM 
framework [63], we have applied it to the partial decay widths of squarks into charginos 
and neutralinos and we have computed the corresponding cross-section. 

In section [2] we present the theoretical framework, by introducing the notation for 
particles and couplings. Section [3J summarizes the results for the one- loop QCD-corrections 
to the squark partial decay widths into charginos and neutralinos. A few details about 
the effective description approach described in [H] are given in this section. This will 
allow readers to know what we have coded into our effective MSSM model for squarks 
within MadGraph. In section HJ we describe our Monte-Carlo tool implementation of the 
MSSM effective couplings. Section [5J presents the numerical analysis: numerical setup 
and parameter choices (section l5.ip . partial decay widths analysis (section l5.2p . and cross- 
section results (section I5T3"|) . Finally section [5] shows our conclusions. 

2 Tree-level relations and parameter definitions 

Here we introduce our notation for SUSY particles and couplings. Throughout this work 
we will use a third-generation notation to describe quarks and squarks, but the analytic 
results and conclusions are completely general, and can be used for quarks-squarks of any 
generation. We will show numerical results only for top-squarks (t), since their signals are 
the most phenomenologically interesting. 

To describe the computation of the partial decay widths, we will follow the conven- 
tions of Ref. [64] . We will study the partial decay widths of sfermions into fermions and 
charginos/neutralinos, 

r(/->/'x). (i) 
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We denote the two sfermion-mass eigenvalues by m? (a = 1,2), with < mj 2 . The 
sfermion-mixing angle Of is defined by the transformation relating the weak- interaction 
(fa = /i> fn) an d the mass eigenstate (f a = fx, $2) sfermion bases: 

f _ d(/) fi . R (f) _ ( c os Of -sinOf \ . , 

~ ' - ^ sin^ cosfl, J ' [Z) 

By this basis transformation, the sfermion mass matrix, 



M| +m} + c 2/3 (T 3 - Q s 2 w ) Mf m f Mf R 

{ ' m f Mf R M? r + m) + Q c 2/3 M | 



becomes diagonal: i?^ A^|-R^^ = diag < m| , m 2 f >. M| is the soft-SUSY-breaking mass 
parameter of the SU(2) L doublets, whereas M 2 - is the soft-SUSY-breaking mass parameter 
of the singlet. T 3 and Q are the usual third component of the isospin and the electric charge 
respectively, mf is the corresponding fermion mass, Mz is the electroweak Z boson mass, 
and Sw is the sinus of the weak mixing angle|§ The mixing parameters in the non-diagonal 
entries read 

M^ R = A b - //tan/3 , M t LR = A t - /i/tan/3 . 

Abj are the trilinear soft-SUSY-breaking couplings, fi is the higgsino mass parameter, 
and tan (3 is the ratio between the vacuum expectation values of the two Higgs doublets, 
tan/3 = V2/V1. The input parameters in the sfermion sector are then: 

(M /L ,M^,M^,A 6 ,A f ,//,tan/3) , (4) 

for each sfermion doublet. From them, we can derive the masses and mixing angles: 

( m b 1 ^ m b 2 ^b) , ( m tv m i2' 6 t) • ( 5 ) 
For the trilinear couplings, we require the approximate (necessary) condition 

A 2 q <3(m 2 + ml + M 2 H + fi 2 ), (6) 

to avoid colour-breaking minima in the MSSM scalar potential [651468] . Here nig is of the 
order of the average squark masses for q — t,b, and M# is the Soft-SUSY-breaking mass 
parameter of the Higgs fields, see [SSHSS] for details. 

Although the tree-level chargino (x + )-neutralino (x°) sector is well known, we give here 
a short description, in order to set our conventions. We start by constructing the following 
set of Weyl spinors: 

r+ = (-iw+,H+) , 

r- = (-iW~,Hr) , (7) 



1 With Mi L = Mi l due to SU(2)l gauge invariance. 

2 We abbreviate trigonometric functions by their initials, like sw = sin 9w, c 2( g = cos(2/3), tw = sw/cw, 
etc. 
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The mass Lagrangian in this basis reads 



C 



M 



-i(r+,r-) 



M T 
M 



r+ 
r- 



-(T 1} T 2} T 3} T 4 )M° 



r 2 
r 3 

vr 4 y 



+ h.c. 



where we have defined 
M = 

M° = 



y/2M w Cp 

M' 
M 

M z cps w -Mzcpcw 

\-M z Sfss w M z s p c w 



Mzcpsw -M z sps w \ 
-M z cpc w M z spc w 





-fi 







(9) 



with M and M' the SU(2)l and U(l)y soft-SUSY-breaking gaugino masses. The four- 
component mass-eigenstate fields are related to the ones in ([7]) by 



Xi 



U*Y~ 



X t 



cx 



u i3 rj 



x c 



cx 



-0T 



where U, V and N are in general complex matrices that diagonalize the mass-matrices 



N*M°N^ 



diag (M l5 M 2 ) (0<M 1 <M 2 ) , 

diag (M°, M 2 °, M 3 °, M 4 °) (0 < M° < M 2 ° < M 3 ° < M 4 °) . 1 j 



Using this notation, the tree-level interaction Lagrangian between fermion-sfermion- 
(chargino or neutralino) reads 



C 



xff 



a=l,2 r 



Xr/«/' 



g f:Xr[Af ar P L + A^ ar P R )f . 



(/) 



Here we have adopted a compact notation, where /' is either / or its SU(2) l partner for 
Xr being a neutralino or a chargino, respectively. Roman characters a,b . . . are reserved for 
sfermion indices and i,j, . . . for chargino indices, Greek indices a, (3, . . . denote neutralinos, 
Roman indices r,s . . . indicate either a chargino or a neutralino. For example, the top- 
squark interactions with charginos are obtained by replacing / — > t, f — > b, Xr — > Xr ■> 
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1,2. The coupling matrices that encode the dynamics are given by 

_ p(*)y* \ p(*)y* 
/>(*) _ \ p(*)rr 

A w aa = 4 (V2A^2iV a4 - Y^t w R^N al ) , 

/i( fe ) _ p( 6 )rr* \ p( b hr* 

: (-v^A^A^ + Y h R t w R®N al ) , (12) 



,(*>) 



with Yl and Y R b the weak hypercharges of the left-handed SU(2)l doublet and right- 
handed singlet fermion, and A t = m t /(y/2Mw sin (5) and = m&/(\/2Mw cos/3) are the 
Yukawa couplings normalized to the SU (2) l gauge coupling constant g. Note the following, 
each coupling is formed by two parts: the gaugino parts (G), formed exclusively by gauge 
couplings, and the higgsino part (H±), which contains factors of the quark masses, each of 
these parts will receive different kind of corrections. First line from equations ( TT2|) can be 
rewritten as: 

A® d = H+ + G. (13) 
The tree-level partial partial decay width reads 

T Tree = fTree^ ^ f , j = _*L_A( m J- , JW* m),) X 

lD7rm^ Ja 1 

J a 

x {m 2 h - M 2 r - mj,) (l4J r | 2 + l^il 2 ) - 4m f/ M r M (a^Ia^S)] ( 14 ) 
with X(x 2 ,y 2 ,z 2 ) = a/[x 2 — (y — z) 2 }[x 2 — (y + z) 2 } 



3 QCD Corrections 

It is known that QCD corrections to the squark partial decay widths into charginos and 
neutralinos can be numerically large, specially in certain regions of the parameter space [12] . 
An effective description of squark/chargino/neutralino couplings, simple to write and to 
introduce in computer codes, was given in (Hj. The complete one-loop corrections to squark 
partial decay widths are already available [I2J.69], but their complicated expressions are 
not suitable for the introduction in Monte-Carlo programs used for experimental analysis. 
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In this article we present the results of the implementation of the effective description in 
MadGraph's MSSM framework. 

In the following we present few details of the effective description approach as given 
in [33]. In this approach, following hints from Higgs-boson physics, an effective Yukawa 
coupling is defined as: 

yeff = = m b(Q) X eff = "V " = ™, t {Q) 

vi fi(l + Am fe )' ' v 2 v 2 (l + Am t ) 

where m q (Q) is the running quark mass and Am, is the finite threshold correction. The 
SUSY-QCD contributions to Am, are 

Am b QCD = ^mgfitanf3 1(m- bi ,m- b2 ,mg), 

Amf QCD = ^ m a-^p Km^m^mg) ( 16 ) 

where I (a, b, c) is the scalar three-point function at zero momentum transfer: 

_ a 2 b 2 In (a 2 /b 2 ) + b 2 c 2 In (b 2 /c 2 ) + a 2 c 2 In (c 2 /a 2 ) 

The effective description of the squark interaction consist in replacing the tree-level quark 
masses in the couplings ffT2l) by the effective Yukawa couplings of eq. (EES]), and use this 
lagrangian to compute the partial decay width, schematically T Yuk ~ e ^ = T Tree {m e q ^) 
(see [33] for details). This expression contains the large one- loop corrections from the finite 
threshold corrections (1151) . but it also contains higher order corrections. At this point one 
can make a computation that combines the higher order effects (which ignore the effects of 
external momenta) and the fixed one-loop (which ignore the higher order effects). At the 
same time, this will allow us to quantify the degree of accuracy obtained by the effective 
description. A Yukawa-improved decay width computation have been defined in [33] and 
they showed that the effective description using just the Yukawa threshold corrections f fT5|) 
is not enough for the squark decay widths description. The one-loop corrections develop 
a term which grows as the gluino mass m g [55], which is absent in the effective Yukawa 
couplings f lT5|) . Therefore, the QCD corrections to squark decay widths produce explicit 
non-decoupling terms of the sort logm^. To understand those terms a renormalization 
group analysis is in order [13]. They constructed an effective theory below the gluino mass 
scale, which contains only squarks, quarks, charginos, neutralinos and gluons in the light 
sector of the theory, and integrate out the gluino contributions. Then, they found out 
the renormalization group equations (RGE) of the gaugino and higgsino couplings, and 
performed the matching with the full MSSM couplings at the gluino mass scale m g . Only 
the logarithmic RGE effects have been considered, neglecting the possible threshold effects 
at the gluino mass scale. Since the effective theory does not contain gluinos, only the 
contributions from the gluon have to be taken into account. We do not present here details 
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of the renormalization group analysis and we restrict to give the results that are relevant 
for our purpose. 

The squark-chargino running coupling constant as a function of the gauge and Higgs 
boson couplings at the renormalization scale Q, up to 0(a s ), are given by [44J: 

4L(Q)-^KW))fl + ^log^ + GA-^log^ . (18) 



7r nig j \ tt rrig 



The logmg term arises as a factor to each higgsino (H±) or gaugino (G) term, eq. ( Fl3l) . 
Note that the expressions for the higgsino and gaugino couplings are different due to the 
different running of the gauge and Higgs-boson couplings between the scales m g and Q. 
The same can be safely extended to the rest of the couplings in f JT2|) . 

Then, the renormalization group running of the coupling constant can be summarized 
as follows: we can use effective gaugino and higgsino couplings given by 

efftris f a s (Q)\W> f a s {Q). Q 
9 n (Q) = 9 7 T -9 1 — log — 



x a s (rrig)J \ n m~ gi 

where X e ^(Q) are the effective Yukawa couplings defined in ( 1T51) . and (3o is the QCD f3- 
function. At this point, a simple expression for the effective description of squark/ chargino/ 
neutralino couplings is given by eqs. (|T5|) . (|T6|) . (|T9|) . Ref. jH] showed that the effective 
description fTT/9]) approximates the improved computation to within 2 — 5% for large enough 
gluino masses (rrig > 1 TeV). The effects of the log-terms are more visible in the gaugino- 
like channels, where the Yukawa couplings play no role, and the bulk of the corrections 
corresponds to the log-terms. In the higgsino-like channels their importance is less appar- 
ent. After introducing these expressions in computer codes a reasonable description for 
squark decays into charginos and neutralinos is accomplished. They can be used in Monte- 
Carlo generators and other computer programs that provide predictions for the LHC and 
the ILC to improve their accuracy, requiring little computational costs. 



4 Monte-Carlo Tool 

The effective couplings of section [3] have been implemented in a standard Monte-Carlo 
tool to allow their use in multiple processes. We have chosen to include them in the 
MadGraph/MadEvent framework |gQHS21EQ|. 

The program MadGraph [60] automatically generates Fortran code to calculate arbi- 
trary tree-level helicity amplitudes. Its later implementation, the MadGraph/MadEvent 
package [HUE2], anows the computation of physical processes at the partonic or hadronic 
level, and includes interfaces to hadronization routines. Recent additions allow also to com- 
pute the partial decay widths of unstable particles (70]. A number of physical models are 
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supplied with the standard version of MadGraph/MadEvent, aside from the SM, we will be 
using the MSSM implementation [63]. For our computation, we use MadGraph/MadEvent 
5.1.4.8 [62], and we make heavy use of its possibility to extend/modify internal physical 
models. MadGraph accepts inputs in the form of SUSY Les Houches Accord (SLHA) [71] 
file format, which allows easy interfacing with programs computing SUSY spectra and 
SUSY particles decay widths. 

We have modified MadGraph's standard MSSM file, changing the squark-quark-neu- 
tralino/chargino vertices by the ones of expression ffl9l) . Using our model file, we are able 
to compute any physical processes involving these vertices, including the leading radiative 
corrections. 

In section l5\2l we compute the partial decay widths of squarks, using the recent addition 
to MadGraph implemented in |70j. These results are checked against an independent 
implementation of the partial decay widths in programs based in the FeynArts/Form- 
Calc/LoopTools [51J[7SJ[7J] (FAFCLT) chain, which were used in Ref. 03]. Our set of 
FAFCLT-based programs prepares a standard SLHA [7T] input file, which contains the 
SUSY inputs, as well as the SUSY spectrum and mixing matrices. This SLHA file is then 
used as input for both, the MadGraph-based and the FAFCLT-based computation. In this 
way we can make meaningful comparisons of the output of both computations. We find 
that both computations agree among them (see section 15.21) and with Ref. [44] . 

In section 15.31 we use our MadGraph-based programs to compute the production cross- 
section of squarks decaying into chargino/neutralino. We generate an SLHA input file with 
our set of FAFCLT-based programs, which we use as input for MadGraph. MadGraph 
needs the full decay width of on-shell particles for the treatment of resonant propagators, 
while our decay programs (FAFCLT- or MadGraph-based) can only compute the fermionic 
decay channels. We use SDECAY 1.1a [74] for the computation of the bosonic decay 
channels partial decay widths. All along the computation a single SLHA file is used for 
all steps, ensuring consistency among the different inputs/outputs. The computing flow is 
the following: 

1. we use the FAFCLT-based programs to write a standard SLHA input file, which 
includes the SUSY inputs as well as the SUSY particle masses and mixing matrices; 

2. we use this file as input for SDECAY to obtain the squark decay tables at the tree- 
level; 

3. then, we use our set of FAFCLT-based programs to replace, in the above file, the 
squark partial decay widths into charginos/neutralinos, and the full decay width, by 
using the effective coupling approach_|; 

4. this SLHA file is used as an input for the MadGraph/MadEvent cross-section com- 
putation, which uses the input section for its internal parameter setup, and the full 
decay width in the resonant particles propagators. 

3 The result would be the same if we used our MadGraph-based programs, since we have checked that 
we obtain the same results. It is a matter of practical convenience. 
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In this way a cross-section at the parton level in the final state is obtained. Mad- 
Graph/MadEvent interfaces may process further this output through parton showers, 
hadronization and detector simulation packages, however this analysis is beyond the scope 
of the present work. 

5 Numerical analysis 
5.1 Numerical set up 

Our programs are able to perform computations for any MSSM parameter space point. 
They admit SLHA [7T] input for easy interaction with other programs/routines. As an 
example of the effects of the new included terms, we will show numerical analysis for fixed 
values of the SUSY parameters, and make plots by changing one parameter at a time. 
On one side, we choose a set of SUSY parameters as described below (see eqs. ( 1201 and 
( 122]) ) and, on the other side, SUSY parameters are chosen from the Snowmass Points and 
Slopes (SPS) [75] set. It should be noted that some of the SPS scenarios are already in 
conflict with LHC data [761177] , but we have chosen them as reference points to perform our 
computation since they have been studied at length in the literature. The results for the 
decay widths serve as a test of correct implementation, they are checked against previous 
computer programs [H]. 

For the SM parameters we use m t = 172 GeV, nib = 4.7 GeV, a s {Mz) = 0.1172, 
s 2 w = 0.221, M z = 91.1875 GeV, 1/at = 137.035989. The renormalization scale Q is taken 
to be the physical mass of the decaying squark. 

First, we take for the central values of the parameters: 

tan/3 = 5 , ll = 300 GeV , M = 200 GeV , M /l = 800 GeV , m g = 3000 GeV , 

M SUSY = M /r = 1000 GeV , A t = A b = 2M /l + fi/ tan (3 = 1660 GeV , ^ ' 

where we have introduced a parameter Msusy as a shortcut for all the SUSY mass param- 
eters which are not explicitly given. We use the grand unification relation M' = 5/3 t^M 
for the bino mass parameter. The value of the trilinear couplings A^t is given by the 
algebraic expression, the given numerical value corresponds to the default values of the 
other parameters, this numerical value will change in the plots, the chosen expression al- 
lows to show plots with a significant parameter variation avoiding colour-breaking-vacuum 
conditions ([6]). The gluino mass is chosen to be large to enhance the effects of the loga- 
rithmic terms. Note also that if the gluino decay channel is open, it will be the dominant 
decay channel for squarks, rendering the chargino/neutralino channels phenomeno logically 
irrelevant. Therefore our region of interest is: 

nig + m q > nig . (21) 
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With these input parameters, the central values for the physical SUSY particle masses are: 

M x+ = (170.40, 337.50) GeV , 

M x o = (89.52, 172.28, 305.46, 338.58) GeV , 

m~ b = (802.05, 1000.30) GeV , 

m~ t = (720.00, 1084.87) GeV . (22) 

Of course, the parameters in eq. (|20|) are just an example for illustrative purposes. We 
denote this parameter choice Ll Def" in the numerical analysis below. In addition we will 
be analyzing the following parameter space points from the SPS [75] collection: lb, 3, 7 
(see below). The corresponding values of the parameters are reproduced in Appendix El 
We have checked that our conclusions hold for a wide range of the parameter space. 

At present the ATLAS collaboration has already put some limits on the top-squark 
mass from pair production ^9TT2Tll24l l25]. the largest excluded top-squark mass at 95% 
C.L. is ~ 465 GeV jT9] assuming BR(ti — > t%°) = 1. Our top-squark mass parameter 
choices are larger than the excluded ones, moreover the exclusion limits would be loosened 
by allowing the existence of several decay channels for the top-squark. 



5.2 Partial Decay Widths 

Following the computation and setup of the previous sections, we perform a complete 
numerical analysis. Here we present numerical results for the partial decay widths of top- 
squarks into charginos/neutralinos by using both the tree-level and one-loop computations 
with FAFCLT-based programs and the MadGraph computations. First, we concentrate 
on the discussion of agreement between the two methods of calculations. Then, we deal 
with the results for the partial decay width in the scenarios defined by the SPS [75] and 
we present a discussion of the accuracy of the effective description of squarks interactions 
in each case. The presentation for the results in the Def scenario, eq. ( 120]) . will be brief, 
since it was largely explored in Ref. [4"4"] . 

First, we present in Figfj] an example of the results of the partial decay width for the 
Def set of parameters (|2"U|) . with the spectrum as described in (|2"2"|) . Here we present the 
results for the two methods of computation, by using MadGraph (labelled MG) and with 
the FAFCLT-based programs (labelled FL), and two approximations: tree-level (labelled 
Tree) and the effective coupling approach of eq. (fl9l) (labelled Eff). Our intention is to show 
explicitly the agreement between the two methods of computation. Specifically, we show 
the partial decay width of the lightest top-squark t\ into charginos and neutralinos as a 
function of tan /3 and m- g . For the purpose of illustration we choose here some of the possible 
decay channels, t\ — > bxt, h — > an d h — > t\2- We focus on these decays because they 
are the ones used to search for top-squarks at the LHC [3| 14], fTT |[T9H2T], [24 1 125]. We see 
perfect agreement between both methods of computation. We have checked explicitly 
that for all possible squarks decay channels into charginos and neutralinos the relative 
deviations between both methods of calculations is within the precision of MadGraph's 
Monte-Carlo numerical integrations (smaller than 0.08%, and typically of the order of 
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Figure 1: Partial decay width of the top-squark t\ into charginos and neutralinos as a func- 
tion of (a) tan/3 and (b) m- g . Shown are the tree-level prediction (Tree) and the effective 
coupling approximation prediction (Eff), in both kinds of computation MadGraph-based 
(MG) and FeynArts/FormCalc-based (FL). Parameters fixed in the Def scenario, eq. ( 120]) . 

0.02%). We have performed a full analysis of the results for the squarks decay widths for 
the Def scenario (120]) and we can conclude that, in the explored range of gluino masses, 
the corrections shrink in more than a 10% the full decay width of squarks compared to 
the tree- level result, and the effect is highly noticed for tan /3 > 35 growing up to a 40% 
in some regions of the parameter space. As expected, the largest difference with respect 
to tree level calculation is achieved at the largest gluino mass evaluated. Results are in 
agreement with Ref. [33]. 

We have made a complete numerical analysis for the Snowmass Points and Slopes 
parameters: SPSlb, SPS3, SPS5 and SPS7. SPS5 provides a very narrow range of useful 
SUSY parameters for our investigation due to a relatively light scalar top quark, therefore 
we are not showing these results in the present paper. The other SPS parameters are 
away of the range of SUSY parameters interesting for our analysis. Notice that the gluino 
mass is around 1 TeV for the three points mentioned above, and therefore the effects of 
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the logarithmic terms is not as enhanced as in the previous analysis, where m g ~ 3TeV. 
We explore the SUSY parameter space for SPSlb, 3 and 7 first by scanning m g in the 
interval [400,5000] GeV while keeping tan/3 fixed, then by scanning tan/3 between 3 and 
50 with rrig fixed. We remind that the best accuracy of the effective description of squarks 
interactions for these points is obtained when m g > 800 GeV. For each set of parameters, 
we compute the partial decay width of top-squarks decaying to charginos and neutralinos 
with both methods of computation and we find a perfect agreement. Figures [2]|3] show a 
summary of the analysis. We analyzed all possible decay channels, however, in Figs. [2j|l]we 
present results for the most interesting channels in the analysis of the cross-section in the 
next section, i.e. the decays of the squarks into the lightest chargino and the two lightest 
neutralinos. We show the tree-level prediction (Tree), the fixed-order one-loop prediction 
(1-loop) and the effective approximation prediction (Eff). In addition to the partial decay 
widths (r) we also show the relative corrections defined as: 

T^correc "pTree/™ \ 

where the label correc is 1-loop for the full one-loop corrections and Eff for the effective 
description computation. 

Figures |2] and [3] show the partial decay widths, T, and relative corrections, 5 ( 123]) . of the 
lightest top-squark t\ decaying into the lightest chargino, (ti — > bxt), and the two lightest 
neutralinos, (ti — > txi 2 ), as a function of m g and tan/3, respectively. In Fig. [2] we see 
that tree-level results are independent of the gluino mass, since the gluino does not enter 
the tree-level process, whereas the effective description follows the logarithmic behaviour 
of the full one-loop corrections. In some figures the value of the relative corrections 5 ( 123]) 
for ti — > txi2 channels are very similar in the SPSlb and SPS3 scenarios, and their lines 
overlap. 

Figure |2] shows that, for the SPSlb and SPS3 scenarios, T{ti — > txt) has positive 
corrections for m g > 850 GeV. The largest difference with respect to tree-level calculation 
is obtained at the largest value of the gluino mass m g = 5TeV (10%). For t\ — > tx^ an d 
ti ~^ bxf decay channels the corrections are negative and slightly decreasing (in absolute 
value) with m g . Opposite to the SPSlb and SPS3 situation, in the SPS7 scenario, X12 
and xf are mostly of higgsino-type and x® a an d xt are mostly of gaugino-type. Then, 
the SPS7 scenario has larger partial decay width values into the lightest chargino and 
neutralinos, and the radiative corrections show a different behaviour. In this scenario all 
decay channels receive negative corrections in the whole analyzed range of gluino mass. 
While the (absolute value) corrections in the t\ — > tx\ channel decrease with m g , the 
t\ — > tx\ and t\ — > bxt channels have an opposite behaviour, becoming more negative as 
m g increases. For all the SUSY parameter choices, the largest partial decay width values 
correspond to the decay channel ti — > bxt- One expects the effective approximation to 
be valid when the gluino mass is much larger than the squark mass scale. In Fig. |2b we 
see that the shape of the effective and one-loop approximation can deviate significantly for 
light gluino masses 400 GeV < m g < 700 GeV, and in this case the effective approximation 
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Figure 2: Lightest top-squark (a) partial decay widths for SPSlb and SPS3 scenarios (first 
row) and SPS7 scenario (second row) and (b) relative corrections, 5 eq. (|23|) . in %, as a 
function of m g . Shown are the tree-level prediction (Tree), the one-loop prediction (1-loop) 
and the effective approximation prediction (Eff). 



is not appropriate. Since the study presented in this work is only relevant when the gluino 
is heavy, we can use the effective approximation since it is valid for m g > 900 GeV. 

Figure [3] shows the same partial decay widths and relative corrections as before, for 
SPSlb, SPS3 and for SPS7, as a function of tan/3. In SPSlb, SPS3 t\ — > tx\ receives 
positive corrections in almost the whole analyzed range of tan/3, with the exception of 
the region of small tan/3 (< 6,13 for SPSlb, SPS3 respectively). For t± — > tx2,bxt 
decay channels, the radiative corrections are always negative. Notice that, while in the 
neutralino decay channels all approximations ( Tree, 1-loop, Eff) to the partial decay widths 
(r) have a similar shape, T(ti — > bxf) behaves differently for each approximation. For all 
SPS at tan /3 w 10, T Tree (ti — > bxt) has a minimum and a steep increase for larger 
tan /3. The radiative corrections are large and negative, compensating this increase. The 
one-loop corrections are proportional to tan/3, and grow faster than the tree-level value, 
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Figure 3: Same as Fig. El but as a function of tan (3. 



the effective computation, on the other hand, resums the large corrections obtained at 
large tan/3, providing a smoother behaviour. For example, in the SPSlb scenario, while 
Y Tree (t l — !. bxt) increases a 25% in the range tan/3 = 20 — 50, the effective prediction 
F E ff(ti — > bxf) increases only a 0.4%. 

For a gluino mass around ~ 900 GeV, corresponding to the nominal value of the chosen 
SPS, the corrections to T{ti — > txi) in SPSlb and SPS3 scenarios are small, below 1% 
in absolute value in the whole tan (3 range. For SPS7 5 is larger (in absolute value), 
evolving from —9% close to —5%. For the x\ decay channel, the deviations 5 are negative: 
-10% < 5 < -7% for SPSlb and SPS3 and 5 w -22% for SPS7. For X i the situation 
changes, 5 l - loov (5 Eff ) reaches -45% (-36%) for SPS7 and -33% (-21%) for SPSlb 
and 3, in absolute values they always increase at a faster rate with tan/3 if compared to 
the neutralino channels. We want to stress that neutralino channels show a nearly flat 
behaviour for 5 because the couplings only involve terms of At strength (oc 1/ tan/3), a 
suppression which wipes out new dependence with tan /3 introduced by Am^ CD and log- 
terms beyond certain tan /3 (^ 10) for the effective approximation. Meanwhile chargino 
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Figure 4: Heavy top-squark (a) partial decay widths for SPS lb and SPS3 scenarios (first 
row) and SPS7 scenario (second row) and (b) relative corrections, 5 eq. (|23|) . in %, as a 
function of rtig. Shown are the tree-level prediction {Tree), the one-loop prediction {1-loop) 
and the effective approximation prediction {Eff). 



channels contain terms of A& strength (oc tan (3) which enhance higher order corrections. As 
pointed above, the SPS7 scenario has the largest radiative corrections. In this scenario, Xi> 
X2 and xt are °f higgsino-type, enhancing higgsino terms (H + ) in the generic couplings of 
expression (flg|) . It is clear that large part of the dependence in tan f3 is encoded in H + , and 
those are also the terms most affected by radiative corrections. SPSlb and SPS3 scenarios 
have the opposite situation, Xii X2 an d Xi are °f gaugino-type, therefore gaugino terms 
{G) -depending much less on tan/3 in the effective approximation- are amplified over H + . 

Figures IU15] show the partial decay widths, T, and relative corrections, 5 eq. fl23|) . for 
the heavy top-squark (t 2 ) decaying into the lightest neutralinos and chargino {Xi2^Xi) as 
a function of the gluino mass (m§) and tan/3, for the parameters SPSlb, SPS3, SPS7. 

Fig. H] shows that the largest partial decay widths correspond to the chargino channel 
(xt) f° r SPSlb/SPS3 and the second neutralino channel (x°) f° r SPS7. The corrections 
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Figure 5: Same as Fig. HI but as a function of tan f3. 



show the characteristic logarithmic behaviour. For SPSlb/SPS3 the corrections are mod- 
erate to large (up to 20%) and they are positive, except for the chargino channel in SPSlb. 
For SPS7 the corrections are negative, they are around —10% for xl, —30% for \2 an d 
—60% for xti m & U cases the effective approximation softens the value of the radiative 
radiative corrections. 

Fig. [5] shows that tan j3 has its largest impact in the chargino xt channel, both in the 
partial decay width and in the radiative corrections. Again, the effective approximation 
moderates the impact of the radiative corrections. Note that, while the one-loop corrections 
surpass the —100% value for the chargino (xt) channel in SPS7 for tan/3 > 41 (meaning 
that the one-loop correction is not reliable), the effective computation predicts a —75%, a 
value that, although quite large, is compatible with perturbation theory, and we can rely 
on this computation. 

Finally, Table [1] summarizes the results for the tree- level computation, Y Tree [GeV] , 
and the relative radiative corrections, 5 E ^ [%] (1231 . of all surveyed partial decay widths 
computed in the effective approximation. Numbers not shown correspond to kinematically 
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0.30750 
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2.1033 
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1.02612 
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4.9 
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2.94559 
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2.34275 
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-54.9 


7.5313 


-16.9 
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4.90028 


-51.7 
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-28.8 
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-16.7 
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Table 1: Tree-level decay width, T Tree [GeV], and relative radiative corrections, 5 E ^ 
[%} (123]) . for the Def (with two values of m g = 1, 3 TeV) and various SPS scenarios. 



closed channels. The results are shown for the three SPS scenarios analyzed in this paper 
SPSlb, SPS3, SPS7, and Def, eq. ( |20l) . We include two sets of Def parameters: one with 
rrig = 3 TeV as in eq. (1201) . and one with m g = ITeV. The reason is a fair comparison 
among the various SUSY scenarios: since the radiative corrections grow as the logarithm 
of rrig, the corrections in the Def parameter set are enhanced over those of the SPS, which 
have a m g ~ 1 TeV. A quick look at Table Q] shows that the radiative corrections are 
quite different in each channel, which means that they will have an impact on the decay 
branching ratios (and production cross-sections). As an example, in SPSlb the tree-level 
prediction is that the largest branching ratio of the heavy top-squark corresponds to the 
second chargino (t 2 — > bxt), followed by \°a an d Xii however, after computing the radiative 
corrections one finds that the largest branching ratio corresponds to Vfo — > txl) — 2-8) 
followed by T(t 2 — > bxt) — 2-6, and T(t 2 — > bxt) — 2-4 in third place. 



5.3 Cross-Section Computation 

After successfully implementing and checking our effective MSSM couplings in MadGraplfl, 
we have computed the cross-section of top-squark pair productions for the Def SPSlb, 
SPS3 and SPS7 scenarios. We have focussed on reactions where top-squarks decay into 
the two lightest neutralinos (x? 2 ) an d the lightest chargino (xf), since these are the search 
channels used at the LHC [3],|1J[TT]. The ATLAS collaboration has already used them to 
analyze LHC data and set limits on top-squark pair-production [T91 - [2T| [2 ^ 123] . 

We have simulated proton-proton (pp) collisions at 14 TeV using the CTEQ6L set of 
parton distribution functions. Our aim is to use the standard MadGraph input files as far as 
possible. By default, MadGraph sets renormalization and factorization scales dynamically. 

4 Files are available on request. 

5 Gauge invariance requires to consider pair production together with single production and non- resonant 
diagrams, see below. 
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(d) (e) (f) 




Figure 6: Generic partonic Feynman diagrams contributing to a{pp — >■ (<zx)(?x))j ( a )-(d) 
double resonant diagrams (cr(pp — >■ g a g* — >■ (<zx)(?x)))> ( e )~(h) single resonant diagrams, 
(i) non-resonant diagram. 

For a production of a pair of heavy particles, as it is our case, these scales are set equal to 
the geometric mean of M 2 + p\ of both particles, where M is the mass of the particle and 
Pt is its transverse momentum. The scale Q of the effective couplings ( I19p is fixed at the 
mass of the internal top-squark for each process. 

In this section we present the results for the top-squark pair production cross-section in 
pp collisions, followed by the decay into a quark and charginos and neutralinos, including 
the effect of the radiative corrections in the effective coupling approximation. 

5.3.1 a(pp -> (q'xr)(q"xs)) 

The cross-section we consider in this work ig^| 

a (pp(-> [q a , q a q* a ]) -> {<1 XrW Xs)) , % = k,2 , q, q" =t,b , Xr,s = Xi, Xi,2 > ( 24 ) 

that is: chargino or neutralino production, associated with top and bottom quarks, which 
proceeds through a single or double top-squark in the intermediate states. Figure [6] shows 

6 Our programs are prepared to compute any final state \r,s — X12 1X1234- 
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the partonic Feynman diagrams. Figs. E^i-d show the double-resonant diagrams contribut- 
ing to the process we are interested in: a(pp — > q a ql — > (qx)(<ix))i however, this set of 
diagrams is not gauge invariant. The complete set of gauge-invariant diagrams includes 
also single- resonant diagrams (Fig. [6fe-h) and non- resonant diagrams (Fig. [6}). We have 
checked that quark- ant i- quark channels (qq, Fig. |6H,h) are three times smaller than the 
gluon-gluon channel (gg), i.e. they represent less than a 30% of the total cross-section. 
Both channels are included in our computation. Since all top-squarks contribute to the 
same final state, one should add up all the amplitudes in Fig. [6] for ti and t 2 , however in 
squark search studies it is costummary to separate the t\ and t 2 channels, and look for each 
signal separately. For this reason we will show results for the cross-section with only one 
of the top-squarks in the intermediate states of Fig. [6j This separation is possible because 
we have checked that the interference effects between t\ and t 2 channels is less than 0.8%. 
We will denote this process as: 

PP -> [%} -> X 

where q a is the squark appearing in the intermediate state, and X is the final state, under 
the understanding that all diagrams (double- resonant, single-resonant, and non-resonant) 
contribute. 

Before performing the full simulation it is useful to analyze approximations to the quan- 
tity under study. Under the narrow width approximation the total cross-section (l24l) can 
be computed as the squark pair production cross-section followed by the decay branching 
ratios, 

a (pp ->■ [q a ] -> {qXrWxslj ^ o(j>p ->- q a q* a ) x BR(q a q' Xr ) x BR(q* a ->- q" Xs ) . (25) 

In this approximation only the double-pole diagrams (Fig. EK-d) contribute. In the present 
process the first part a{pp — > qiq*) proceeds purely through standard QCD couplings, the 
only SUSY parameters being the squark masses. All the information on SUSY couplings is 
in the second part of the expression, that is, in the decay branching ratios or partial decay 
widths, which have been analyzed in section 15.21 The top-squark production cross-section 
has been computed to NLO-SUSY-QCD [78] and NLL-SUSY-QCD [79], the radiative cor- 
rections increase the production cross-section at the 14 TeV LHC by a 25%— 50% depending 
on the top-squark mass, with a mild dependence on other SUSY parameters [781179]. 

In order to asses the effects of the radiative corrections, it is useful to define a factor 
Ksusy as: 

7^ ^Effective fr .„^ 

J^SUSY — , (/DJ 

CTree 

which under the narrow-width approximation of eq. ( 125|) can be computed as a ratio of 
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partial decay widths: 



K 



SUSY 



Kpartiai 
Kf u ll 



partial — ^ ~ : ~ X 



T 1 , 



Y E ff 

K fu ii = I I ( 27 ) 



where is the full decay width and Eff and Tree stand for the effective approximation 
and the tree-level computation. For the parameter space points of the present study, all 
bosonic t\ decay channels are kinematically closed, and therefore the full decay width 
consists only of the chargino/neutralino (inos) channels. Bosonic channels have some 
impact on the ti full decay width. Ksusy is composed of two factors: K part iai, regarding the 
partial decay widths of the produced squarks decaying into a quark and the selected final 
particles (x? 2 or Xi analyzed in section 15.2 j) and Kf u u the radiative correction to the total 
decay with. For the parameter space explored in the present work we have found that the 
radiative corrections to the total decay width are always negative (Kf u u < 1). Of course, 
using the narrow-width approximation (I25p one can not compute angular distributions or 
correlations, something that it is possible to perform with the MadGraph implementation. 

The results for the cross-sections are summarized in Figures I71fl~2| where we show the 
tree-level and effective approximation cross-section computations, as well as the Ksusy (J2"6]) 
factor, function of tan/3 and m~ g . 

Figures HE] present the cross-section corresponding to the lightest top-squark t±. Fig.[7K 
shows the results for a{pp — > [tx] — > (tXx)(tXx)) as a function of rrig and tan/3 for the 
different scenarios of the SUSY input parameters analyzed in this work. The ratio Ksusy 
is shown in Fig. [7b. The radiative corrections are positive in all scenarios, this is due 
to the combination of two factors: first, the radiative corrections to the partial decay 
width r(?i -> 6x?) 

are positive in the scenarios SPSlb, SPS3 and Def (Figs. [2, [3] and 
Tabled]), providing a K partia i > 1, but, more importantly, the radiative corrections to the 
total decay width are negative, providing a Kf uit < 1, and therefore a K S usy > 1 (see 
eq. (j2~Tj) ). as can be seen in Table [TJ For example, in SPSlb the leading decay channels (at 
the tree-level) are bxf and b\t > but they have negative radiative corrections of —13.7% 
and —39.9% respectively, in the end, when taking into account all channels one obtains 
Tf u f = 4.49 GeV, V E f H = 3.78 GeV, providing K fuil = (0.84) 2 ~ 0.71 ~ 1/1.41, so a 
fifuii = ~ 26% radiative correction to the total decay width translates to a 41% enhancement 
of the production cross-section. The increase of the radiative corrections Ksusy with m- g 
is given by the log(mg)-terms ffl9l) of the effective computation. This figure shows the 
importance of the radiative corrections, which can enhance the cross-section by a large 
factor (between 1.2 and 5 in Fig. [7J, and also of the newly included log(mg)-terms, the log 
terms produce an increase of the KsusY-f&ctor (and hence the cross-section) of 1.6/1.4 = 
1.14, 1.4/1.2 = 1.17, 2/1.5 = 1.33, 2.4/1.7 = 1.41 in the range m- g = 1000 - 5000 GeV for 
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Figure 7: (a) Cross-section of pp — > [ti] — >■ tx^tXi an d (b) K S usy, as a function of tan/3 
and rrig for the input parameters as denned in SPSlb, 3, 7 and Def ( 120]) . 



the scenarios SPSlb, SPS3, SPS7, Def respectively. The largest cross-sections are obtained 
in the SPSlb and SPS3 scenarios, which have a similar behaviour of the corrections. At 
tan/3 = 30, the cross-section values in the effective approximation are 4.08 x 10~ 3 pb 
(SPSlb) and 5.38 x 10~ 3 pb (SPS3), a factor 1.43 and 1.37 larger than the tree-level 
prediction respectively. The SPS7 scenario provides an intermediate value of the cross- 
section in our analysis. For tan /3 = 15 and rrig = 926 GeV (as fixed in this scenario) the 
effective prediction for the cross-section is 1.82 x 10 -4 pb, a factor 1.56 larger than the tree- 
level prediction (see Fig. [7b). The tree- level cross-sections decrease as a function of tan/3 
in the region tan (3 > 10, this behaviour is softened by the positive radiative corrections. 
We recall that the squark masses (and hence phase-space factors) also change with tan /3, 
i.e. mi changes from 772 GeV to 780 GeV. The smaller values of the cross-section are 
obtained in the Def scenario (1201) . e.g. a E ^ = 6.54 x 10~ 6 pb for the effective computation 
prediction at tan/3 = 5, a factor 2.25 larger than the tree-level prediction. The largest 
radiative corrections are obtained at large tan/3, with Ksusy = 4.81. We recall that in 
this scenario rrig = 3 TeV, enhancing the radiative corrections. 

The results for a{pp — > [ij — > tx^tXi) are presented in Fig. [HI For SPSlb/SPS3 
the cross-sections are slightly larger than the (tXx)(^Xi) channel (Fig. [7]), whereas they 
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Figure 8: (a) Cross-section of pp — > [ti] — >■ tx^pXi an d O 3 ) K S usy, as a function of tan/3 
and for the input parameters as denned in SPSlb, 3, 7 and Def ( 120]) . 



are much larger for SPS7 / Def. We can trace back this behaviour to the relative value 
of the partial decay widths (Figs. El |3l Table [[)) and combinatorial factors. First of 
all, there is a factor 2 combinatorial enhancement factor because now the particles in 
the final state are different. Then in SPSlb/SPS3 there is a slight suppression because 
BR(ti — > tx%) < BR{ti — > txi), whereas in SPS7/ Def there is a large enhancement be- 
cause BR(ti — > tx%) ^ BR(fi — > txi). The radiative corrections (Fig. |Hb) are smaller than 
in the (tXi)(*Xi) channel (Fig. [7b), this is due to the fact that the radiative corrections 
are smaller (more negative) for the partial decay width V(t\ — > tx®) than for T(ti — > tx\) 
(Table [1]). This channel in the Def scenario has a cross-section one order of magnitude 
larger than the (tXiX^Xi) channel. 

To finish with the t\ production analysis, Fig. [H] shows the cross-sections and radiative 
corrections for the lightest neutralino (xl) an d chargino (xf) production mediated by 
t\. Since the experimental analysis does not perform charge identification, both charge- 
conjugate channels a{pp — > [ti] — > tx^bXi) an d o{pp — > [ti] —> ixibxt) contribute to the 
same experimental signal, therefore Fig. |HK shows the sum of both cross- sect ions^. This 

7 Bottom-squark mediated channels also contribute to this final state, we don't take into account their 
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Figure 9: (a) Cross-section a{pp — > [ti] — > bxftXi) + cr(pp —> \px] — > b\i tXi) and (b) 
Ksusy, as a function of tan/3 and trig for the input parameters as defined in SPSlb, 3, 7 
and Def flU. 



channel provides the largest cross-section for ^-pair production in all studied scenarios. 
For example, at tan/3 = 30, the cross-section values in the effective approximation are 
1.89 x 1(T 2 pb (SPSlb), 2.29 x 10~ 2 pb (SPS3), 1.85 x 1(T 3 pb (SPS7) and 3.08 x 10~ 4 
pb (Def). The radiative corrections (Fig. [9b) are slightly smaller in this channel for the 
SPS scenarios, providing Ksusy factors in the range Ksusy — 11 — 1-4 (10% — 40% 
increase of the cross-section). The Def scenario has the largest radiative corrections for 
the reference value tan j3 = 5 f!20p , but they have a different behaviour than the previous 
channels (Figs. El EJ) at large tan/3: they decrease instead of growing at large tan/3. In 
this scenario the leading t\ decay channels are X34 and xti & U of them having negative 
corrections (Table [1]) which grow with tan /3. This provides a growing contribution to 
Ksusy through 1/Kf u u (T27j). The other contribution, K partia i, has two distinct factors: 
the radiative corrections to T(ti — > tXi) have a nearly flat behaviour with tan at a value 
around 18% — 22%, but the radiative corrections to F(ti —> bxf) have a strong decreasing 

contribution because they have different kincmatical resonances, and the interferences with t channels are 
small. 
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Figure 10: (a) Cross-section of pp — > [t 2 ] — > tX\tx\ and (b) K S usy, as a function of tan/3 
and rrig for the input parameters as denned in SPSlb, 3, 7 and Def ([20]) . 

behaviour from ~ 29% (tan/3 = 5) to ~ —52% (tan/3 = 50). This decrease in K partia i 
partially compensates for the increase due to 1/Kf u u resulting in a more flat behaviour 
with tan /3 as compared with the tXatX°p channels. 

The results for the t 2 mediated cross-sections are shown in figures [T0lfT2l The values of 
these cross-sections are smaller than for the t\ mediated cross-sections. Fig. [10] presents the 
cross-section a(pp — > \t 2 } — > txlixi) and the corresponding radiative corrections. Radiative 
corrections are positive (except for a small corner at low tan/3 for SPS7), enhancing the 
cross-section, which stays in the 10~ 7 — 10~ 5 pb range. The largest results for the cross- 
section are obtained in the SPSlb and SPS3 scenario, which have similar values, and they 
overlap in the plots as function of tan (3. At tan f3 = 30, the corresponding values in the 
effective approximation are 1.094 x 10~ 5 pb (SPSlb) and 1.098 x 10" 5 pb (SPS3). The 
different results between these points (when plotted as a function of rrig) arise from the 
different nominal values of tan/3 (appendix [A}. The same similarities occur in the other 
channels (Figs. [TT] [T2"]) . but not as pronounced as the present one. In scenario Def we 
obtain an intermediate value of the cross-section (3.37 x 10 -6 pb), and the lowest value is 
obtained in the SPS7 scenario (7.29 x 10 -7 pb). Note the quite distinct behaviour of the 
radiative corrections (Fig. [TUb) as compared with the ti channels (Fig. [7b): in the present 
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Figure 11: (a) Cross-section of pp — > [t 2 ] — >■ tx^Xi an d O 3 ) K S usy, as a function of tan/3 
the input parameters as denned in SPSlb, 3, 7 and Def (!20|) . 



and for 



channel the Def scenario has smaller radiative corrections, and a flatter evolution with 
tan/3. Quite opposite, the SPS scenarios have larger radiative corrections and a steeper 
slope as a function of tan/3. The radiative corrections tend to soften the slopes of the 
cross-sections as a function of tan /3. We note also a small region of negative radiative 
corrections for SPS7 at low tan < 7. 

Figure [UJ shows the results for a(pp — > [£2] — > tx^tXi)- The cross-sections are larger 
than in the X1X1 channel in all scenarios, being in the range of 10~ 6 — 10 -5 pb. The largest 
values are again obtained for the SPSlb and SPS3 scenarios. In these two scenarios, the 
partial decay widths of I2 decaying into \\ 2 have always positive radiative corrections in the 
whole explored region of gluino masses and tan /3 (Figs. Hb, [5b), providing an enhancement 
to Ksusy-, additional to the negative correction to the total decay width [Kf u u < 1). 
The cross-sections for SPS7 and Def have a mild evolution with tan/3, and the radiative 
corrections are smaller than in the X1X1 channel. 

Finally, we present the results for the chargino-neutralino channel in Fig. [12j Again, 
the plots in Fig. [12] contain the sum of the two charge-conjugate modes: a{pp — > [t 2 ] -> 
bxttXi) + &(pp — > [t 2 ] — > bxitxi)- This channel provides the largest cross-section for 
t 2 production in the SPSlb, SPS3 and Def scenarios. In these scenarios the radiative 
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Figure 12: (a) Cross-section a{pp — > [t 2 ] —> bxt^Xi) + &{pv ~^ IP2) ~^ &Xi ^Xi) an d (b) 
Ksusy, as a function of tan/3 and m g for the input parameters as defined in SPSlb, 3, 7 
and Def fl20]). 



corrections are smaller than previous channels, and they have a quite flat evolution with 
tan/3 and rrtg. In contrast, in the SPS7 scenario we find a steep evolution with tan/3 of 
both: the cross-section and the radiative corrections. This evolution is inherited from the 
partial decay width T(t2 — > bxt) (Fig. [5]). The radiative corrections are negative in a wide 
region of parameters, producing a decrease of the cross-section. 

Let us note that the main t 2 decay channels correspond to X3 4 m the neutralino sector 
and xt i n t ne chargino sector (Tabled]), together these channels are responsible for more 
than 60% of the t 2 decay branching ratio in all studied scenarios. At the same time, these 
leading channels have large negative radiative corrections (l^^l > 20%), providing a large 
1/Kfuu f )27|) . For the SPSlb, SPS3, Def scenarios the radiative corrections to T(t2 — > bxt) 
are quite moderate and negative (Fig. [51 Table [[]) producing a small reduction of Ksusy 
through K partia i. But for SPS7 the radiative corrections are larger, and have a very strong 
tan/3 dependence 40% > 5 Eff > -75% for 5 < tan/3 < 50 (Fig. ED, for tan/3 > 10 the 
negative corrections to T(t 2 —> bxt) (negative contribution to Ksusy) overcompensate 
the negative corrections to the total decay width Tf u u (positive contribution to Ksusy), 
providing an overall negative contribution to the cross-section radiative corrections. 
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Table 2: Effects of the radiative corrections to production cross-sections in the effective 
approximation, K SUSY eq. (|26|) . for SPSlb, SPS3, SPS7 and Def SUSY parameters choice 
at different tan /3. 

To complete the discussion in this section, Table |2] shows the radiative correction factor 
Ksusy Q26p for all the production channels cross-sections and all the scenarios presented 
in this work. Results are shown for three different values of tan /3, to make the comparison 
more meaningful. For the same reason the Def scenario is shown also for a light gluino mass 
m g = 1 TeV. In general, the radiative corrections are positive, increasing the production 
cross-section. As analyzed above, this is mainly (but not only) due to the fact that the 
radiative corrections decrease the top-squark total decay width. For m g = 1 TeV the 
positive corrections provide a K S usy factor between 1.05 and 3.91, depending on the 
channel and scenario. Negative radiative corrections are only obtained for a(pp — > [t 2 ] — > 
(bxt)(tXi)) i n the SPS7 scenario and the Def scenario at m g = 1 TeV, due to large negative 
corrections to r(t 2 — > bxt) (Fig. [5] and Tabled]), they provide a Ksusy between 0.62 and 
0.99 depending on the scenario. The largest radiative corrections are obtained at large 
tan/3. We recall that the corrections grow with the gluino mass (e.g. Ksusy = 4.81 for 
tan/3 = 50, m-g = 3 TeV in the (xiXi) channel and Def scenario). 

6 Conclusions 

We have implemented and tested an effective description of squark interactions with 
charginos and neutralinos in the MSSM [44] into MadGraph [60H62] . A careful check 
of our implementation has been done by comparing the computation of the partial decay 
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width of squarks into charginos and neutralinos with FeynArts/FormCalc/Loop Tools-based 
programs [6lJ[72j[73]. We find perfect agreement. We have reproduced previous results in 
the literature [44]. This implementation allows to perform any Monte-Carlo computation 
taking into account the leading SUSY radiative corrections to squark-chargino/neutralino 
couplings. 

Using this implementation, we have computed the partial decay widths of top-squarks 
into charginos and neutralinos in the effective coupling approximation, for sets of SUSY 
input parameters as defined in the Snowmass Points and Slopes SPSlb, SPS3 and SPS7, 
which correspond to scenarios in which the effective approximation can be applied, eq. (1211) . 
We have also analyzed a particular scenario which we denote as Def, eq. ( 1201) . We have 
compared the effective approximation to the fixed-order one-loop computation, and find 
that the effective approximation provides a good description of the radiative corrections if 
the gluino mass is larger than the squark mass (jrig > 800 GeV for the chosen scenarios). 

For the first time we have computed SUSY particle production cross-sections at the 
14 TeV LHC using the effective description of squark interactions. We have analyzed the 
cross-sections performing a comprehensive scan of the SUSY parameter space around the 
scenarios SPSlb, SPS3, SPS7 and Def We have focussed on reactions involving top- 
squarks, giving rise to a final state involving the lightests neutralinos and chargino (x?2) 
Xi)- These are the channels used by the CMS and ATLAS collaborations to perform 
squark searches. The radiative corrections are positive in most of the explored parameter 
space, producing an increase of the SUSY production cross-section. For a gluino mass of 
nig ~ 1 TeV they provide up to a factor 4 enhancement of the cross-section (Table [2]). This 
enhancement factors are mostly (but not only) driven by the negative radiative corrections 
to the top-squark total decay width Kj uU (1271) . The corrections grow with the gluino mass 
(m g ). This leads to the lucky situation that, if the gluino is heavy (and hence, has a small 
production rate at the LHC) the radiative corrections to the squark-chargino/neutralino 
couplings will be large, and easier to study at the LHC. On the other hand, if the gluino 
is light the radiative corrections to the squark-chargino/neutralino couplings will be small, 
but then the gluino production rate at the LHC will be quite large, and will be easy to 
study! all in all, it is a win- win situation (provided SUSY exists at all). We leave the prove 
or rejection of SUSY to our experimental colleages, and give them a new tool to explore 
with better precision the SUSY parameter space at the LHC or the future ILC. 
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A Snowmass points and slopes parameters 

Snowmass ponts and slopes paremeters from Ref. [75J. Mass parameters in GeV. 
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